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This paper is about the stagnation point flow and mass transfer with chemical reaction past a stretching/ 
shrinking cylinder. The governing partial differential equations in cylindrical form are transformed into 
ordinary differential equations by a similarity transformation. The transformed equations are solved 
numerically using a shooting method. Results for the skin friction coefficient, Schmidt number, velocity 
profiles as well as concentration profiles are presented for different values of the governing parameters. 
Effects of the curvature parameter, stretching/shrinking parameter and Schmidt number on the flow and 
mass transfer characteristics are examined. The study indicates that dual solutions exist for the shrinking 
cylinder but for the stretching cylinder, the solution is unique. It is observed that the surface shear stress and 
the mass transfer rate at the surface increase as the curvature parameter increases. 

The stagnation flow is about the fluid motion near the stagnation point. The fluid pressure, the heat transfer 
and the rate of mass deposition are highest in the stagnation area. Wang 1 investigated the stagnation flow 
towards a shrinking sheet and found that the convective heat transfer decreases with the shrinking rate due 
to an increase in the boundary layer thickness. He obtained dual solutions and unique solution for specific values 
of the ratio of shrinking and straining rates. The similar flow over a shrinking sheet in a micropolar fluid was 
investigated by Ishak et al. 2 . Bachok et al. 3 considered the stagnation-point flow and heat transfer from a warm, 
laminar liquid flow to a melting stretching/ shrinking sheet. This problem was then extended to a micropolar fluid 
by Yacob et al. 4 . The effects of homogeneous-heterogeneous reactions on the steady boundary layer flow near the 
stagnation-point on a stretching/shrinking surface were studied by Bachok et al. 5 . Later, Bachok et al. 6,7 studied 
the boundary layer stagnation-point flow towards a stretching/shrinking sheet in a nanofluid. Bhattacharyya 8 
discussed the existence of dual solutions in the boundary layer flow and mass transfer with chemical reaction 
passing through a stretching/shrinking sheet. He found that the concentration boundary layer thickness decreases 
with increasing values of the Schmidt number and the reaction-rate parameter for both solutions. The problem 
considered by Wang 1 was extended by Bhattacharyya 9 , Bhattacharyya and Layek 10 , Bhattacharyya and Pop 11 , 
Bhattacharyya et al. 12,13 and Lok et al. 14 to various physical conditions. 

The flow over a cylinder is considered to be two-dimensional if the body radius is large compared to the 
boundary layer thickness. On the other hand, for a thin or slender cylinder, the radius of the cylinder may be of the 
same order as that of the boundary layer thickness. Therefore, the flow may be considered as axi-symmetric 
instead of two-dimensional 15,16 . The study of steady flow in a viscous and incompressible fluid outside a vertical 
cylinder has been done by Ishak 17 and Bachok and Ishak 18 . The effect of slot suction/injection as studied by Datta 
et al. 15 and Kumari and Nath 16 may be useful in the cooling of nuclear reactors during emergency shutdown, 
where a part of the surface can be cooled by injection of coolant (Ishak et al. 19 ). Lin and Shih 20,21 considered the 
laminar boundary layer and heat transfer along horizontally and vertically moving cylinders with constant 
velocity and found that the similarity solutions could not be obtained due to the curvature effect of the cylinder. 
Ishak and Nazar 22 showed that the similarity solutions may be obtained by assuming that the cylinder is stretched 
with linear velocity in the axial direction and noted that their study is the extension of the papers by Grubka and 
Bobba 23 and Ali 24 , from a stretching sheet to a stretching cylinder. 

The addition of chemical reaction in the boundary layer flow has huge applications in air and water pollutions, 
fibrous insulation, atmospheric flows and many other chemical engineering problems. Hayat et al. 25 discussed the 
mass transfer in the steady two-dimensional MHD boundary layer flow of an upper-convected Maxwell fluid past 
a porous shrinking sheet in the presence of chemical reaction and expressions for the velocity and concentration 



SCIENTIFIC REPORTS | 4:4178 | DOI: 1 0. 1 038/srep041 78 



1 



profiles were obtained using HAM. Bhattacharyya and Layek 26 dis- 
cussed the behavior of chemically reactive solute distribution in 
MHD boundary layer flow over a permeable stretching sheet, vertical 
stretching sheet 27 and stagnation-point flow over a stretching sheet 28 . 
The aim of the present study is to extend the paper by Bhattacharyya 8 
to a cylindrical case. We investigate the skin friction and the mass 
transfer characteristics at the solid-fluid interface in the presence of 
chemical reaction. To the best of our knowledge, this problem has not 
been studied before and all results are new. 



where r is the coordinate measured in the radial direction, and u and 
v are the velocity components in the x and r directions, respectively. 
Further, T is the temperature in the boundary layer, v is the kinematic 
viscosity coefficient, C is the concentration, C m is the constant 
concentration in the free stream, D is the diffusion coefficient and 
_R r denotes the reaction rate of the solute. 
The boundary conditions are 

u = u w = cx/L, v = 0, C = C W at r = R, 



Problem formulation. Consider a steady stagnation-point flow 
towards a horizontal linearly stretching/shrinking cylinder with 
radius R placed in an incompressible viscous fluid of constant 
temperature T w and chemically reactive species undergoing first 
order chemical reaction as shown in Fig. 1. It is assumed that the 
free stream and the stretching/shrinking velocities are u e = ax/L and 
u w = c x/L respectively, where a and c are constants, x is the 
coordinate measured along the cylinder and L is the characteristics 
length. The boundary layer equations are (Ishak 17 ; Bhattacharyya 8 ) 



d d 
1 r(ru)+ — (rv) = 0, 
ox or 



8x 



du 

-v—=u e 
or 



dx 



+ v 



dC 8C 
ox or 



D 



8 2 C 
'dr 2 ' 



d 2 u 1 du 
dr 2 r dr 



18C 
r dr 



-R r (C-C a 



(1) 



(2) 



(3) 



>u e = ax/L, C— >Cqo as r- 



(4) 



We look for similarity solutions of equations (l)-(3), subject to the 
boundary conditions (4), by writing 



f i-R 2 
2R 



1/2 /va\ 1/2 r/ ^ C-C a 



C w Co; 



where ;] is the similarity variable, \jj is the stream function defined as 
u — r^ 1 dip /dr and v= — r~ 1 dtft/dx, which identically satisfies equa- 
tion (1). By defining r\ in this form, the boundary conditions at r = R 
reduce to the boundary conditions at v\ = 0, which is more conveni- 
ent for numerical computations. 

Substituting (5) into equations (2) and (3), we obtain the following 
nonlinear ordinary differential equations: 



(l+2y )? )/'" + 2 > /"+#" + l-/' 2 =0, 



(6) 




Figure 1 | Physical model and coordinate system for cylinder (a) stretching case (b) shrinking case. 
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[ ] second solution. 



( 1 + 2y>,)f' + 2>f + Scf(t>' - Scfi<l> = 0, (7) 
subject to the boundary conditions (4) which become 

/(o)=o, /'(o)=s, <Ko)=h 

/'(oo)^l, ^(oo)^0, 



(8) 



where y is the curvature parameter, Sc is the Schmidt number, fS is the 
reaction-rate parameter defined respectively as 



vL 



1/2 



D a 



(9) 



and e = c/a is the stretching/shrinking parameter with e > 0 is for 
stretching and e < 0 is for shrinking. 

The main physical quantities of interest are the value of/"(0), 
being a measure of the skin friction, and the concentration gradient 



/"(0) 



y=0, 0.2,0.4 




— <l>'(0). Our main aim is to find how the values of/"(0) and —<j)'(0) 
vary in terms of parameters y, Sc and /J. When y = 0 (flat plate), the 
present problem reduces to those considered by Bhattacharyya 8 . 

Results and Discussion 

Numerical solutions to the ordinary differential equations (6) and (7) 
with the boundary conditions (8) form a two-point boundary value 
problem (BVP) and are solved using a shooting method, by convert- 
ing them into an initial value problem (I VP) . This method is very well 
described in the recent papers by Bhattacharyya 8 , Bhattacharyya et 
al. 12 and Bachok et al. 29 . In this method, we choose suitable finite 
values of r\, say which depend on the values of the parameters 
considered. First, the system of equations (6) and (7) is reduced to a 
first-order system (by introducing new variables) as follows: 

f'=p, p' = q, (l + 2y»;)g' + 2yq+fq+ l-p 2 = 0, (10) 
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Figure 2 | Skin friction coefficient f"(0) with e for various values of 7. 



Figure 3 | Variation of — (j)' (0) with £ for various values of y when Sc = 1 
and P= I. 
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Table 2 | Variations of s c with curvature parameter y 

y Bhattacharyya 8 Present study 



0 
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-1.24657 



-1.24657 
-1.38090 
-1.49353 



$ = r, ( 1 + 2yfi)r' + 2yr + Scfr - Scptj) = 0, 
with the boundary conditions 

/(0) = 0, p(0) = 8, <K0)=1, 

Now we have a set of 'partial' initial conditions 

/(0)=0, p(0) = e, q(0) = ?, 0(0) = 1, r(0) = ?. 



(11) 



(12) 



(13) 



As we notice, we do not have the values of q(0) and r(0). To solve Eqs. 
(10) and (11) as an IVP, we need the values of q(0) and r(0), i.e.,/'(0) 
and ^'(0). We guess these values and apply the Runge-Kutta- 
Fehlberg method in maple software, then see if this guess matches 
the boundary conditions at the very end. Varying the initial slopes 
gives rise to a set of profiles which suggest the trajectory of a projectile 
'shot' from the initial point. That initial slope is sought which results 
in the trajectory 'hitting' the target, that is, the final value (Bailey 
etal. 30 ). 

To determine either the solution obtained is valid or not, it is 
necessary to check the velocity and the concentration profiles. The 
correct profiles must satisfy the boundary conditions at r\ = r\ x = 30 
asymptotically. This procedure is repeated for other guessing values 
of q(Q) and r(0) for the same values of parameters. If a different 
solution is obtained and the profiles satisfy the far field boundary 
conditions asymptotically but with different boundary layer thick- 
ness, then the solution is also a solution to the boundary-value prob- 
lem (second solution). This method has been successfully used by the 
present authors to solve various problems related to the boundary 
layer flow (see Najib et al. 31 and Bachok et al. 3,7 ). 
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Figure 5 | Variation of — </>' (0) with e for various values of P when Sc = 0.5 
and y — 0.2. 

Table 1 shows the variations and comparison of/'(0) with those of 
previous researcher (Bhattacharyya 8 ) for the flat-plate case, which 
show a good agreement, thus give confidence that the numerical 
results obtained are accurate. Other than that, the values of/'(0) 
for y = 0.2 and y = 0.4 are also included in Table 1 for future 
references. 

Variations of the skin friction coefficient /'(0) and the concentra- 
tion gradient —<f>'(0) with e and y are shown in Figs. 2 and 3 when Sc 
= 1 and P = 1. From the figures, we can see that there are regions of 
unique solutions for n a — 1.0, dual solutions for 8 C < e < — 1.0 and 
no solutions for s < 8 C < —1.0. Therefore, the solutions exist up to 
the critical value e = e c (< — 1.0). The boundary layer approximation 
breaks down at e = s c ; thus we are unable to obtain further results for 
s < e, c . Beyond this value, the boundary later has separated from the 
surface. The critical value of £ (say £ c ) are presented in Table 2, which 
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Figure 7 | Concentration profiles for various values of /when p=l, Figure „ | Concentration proflles m for various values of p when e , 
£ = —1.2 and Sc — 0.5. 



show a very good agreement with that of Bhattacharyya 8 , for the flat 
plate case (y = 0). It is evident from figures 2-5 and Table 2 that |e c | 
increases with an increase in the curvature parameter y. The range of 
e for which the solution exists is larger for y > 0 (cylinder) compared 
to y = 0 (flat plate). Thus, this demonstrates that a cylinder increases 
the range of existence of the similarity solutions to the equations (6)- 
(8) compared to a flat plate, i.e. the boundary layer separation is 
delayed for a cylinder. The results shown in Fig. 2 also indicate that 
as the curvature parameter y increases, the skin friction coefficient 
/'(0) also increases. Fig. 3 shows the values of concentration gradient 
— <j)'(0) which are proportional to the rate of mass transfer with an 
increase of"/ when Sc = 1 and /? = 1. Moreover, Fig. 4 presents the 
variation of the concentration gradient —(f>'(0) for various values of 
Schmidt number Sc when fS = 0 and y = 0.2. It is seen that when 
Schmidt number Sc increases, the concentration gradient also 



-1.2,7= 0.2 and Sc= 0.5. 

increases. Fig. 5 shows the variation of concentration gradient 
— ^'(0) for increasing fS when Sc = 0.5 and y = 0.2. In general, 
increasing \] is to increase the concentration gradient at the surface. 
The mass transfer increases with increasing values of e, for the first 
solution, but decreases with increasing e for the second solution. 

The variations of velocity and concentration profiles for different 
values of y, £, Sc and ft are display in Figs. 6-11, which show that the 
far-field boundary conditions are satisfied, and thus support the 
validity of the numerical results obtained. The dual velocity profiles 
/ {if) in Fig. 6 show that the velocity increases with increasing y for the 
first solution and conversely for the second solution, it decreases. It is 
to be noted that the momentum boundary layer thickness for the 
second solution is thicker than the thickness of the first solution. On 
the other hand, it should be mentioned that the first solutions are 
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Figure 11 | Concentration profiles (j>{tf) for various values of e when /?= 1, 
y — 0.2 and Sc — 1. 

stable and physically realizable, while the second solutions are not. 
The procedure for showing this has been described by Weidman 
et al. 32 , Merkin 33 , and recently by Postelnicu and Pop 34 , so that we 
will not repeat it here. From Fig. 7, we noticed that the value of 
concentration profile <j)(i]) initially decreases with y and after that 
for large r\, changing the nature it increases with y. In Fig. 8, the effects 
of Schmidt number Sc on the concentration profile <f>(r]) is exhibited. 
The dual concentration profiles of Fig. 8 demonstrate that concen- 
tration decreases with Sc. As a result, the concentration boundary 
layer thickness reduces with enhancement of Sc. Fig. 9 presents the 
concentration profile <f>(r\) at a fixed value of Sc, e, and y for both 
solutions, which shows that the boundary layer thickness decreases 
with increasing p. 

Conclusion 

We have theoretically investigated how the governing parameters, 
i.e. curvature parameter y, stretching/shrinking parameter e and 
Schmidt number Sc, influence the boundary layer flow and mass 
transfer with chemical reaction characteristics on a stretching/ 
shrinking cylinder. It was found that the solutions for the stretching 
case are unique but for the shrinking case, there are dual solutions for 
a certain range of the stretching/ shrinking parameter. The curvature 
parameter y increases the range of existence of the similarity solu- 
tions, which in turn delays the boundary layer breakdown. Hence, 
the boundary layer separation is delayed for a cylinder (y > 0) com- 
pared to a flat plate (y = 0). Further, it was found that increasing the 
curvature parameter y is to increase the surface shear stress and the 
mass transfer at the surface. The concentration boundary layer thick- 
ness decreases with increasing values of Schmidt number and reac- 
tion-rate parameter for both solutions. 
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